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We present a formalism for studying the behaviour of quantum systems coupled to nonequilibrium 
environments exhibiting nonGaussian fluctuations. We discuss the role of a qubit as a detector of 
the statistics of environmental fluctuations, as well as nonMarkovian effects in both weak and strong 
coupling limits. We also discuss the differences between the influences of classical and quantum en- 
vironments. As examples of the application of this formalism we study the dephasing and relaxation 
of a charge qubit coupled to nonequillibrium electron transport through single and double quantum 
dots. 
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The standard paradigm of system-environment inter- 
actions in quantum mechanics employs an equilibrium 
environment which is large enough that its fluctuations 
are Gaussian This model is inappropriate, how- 

ever, if our quantum system couples strongly to a small 
number of environmental degrees of freedom, which will 
typically be out of equilibrium and display a full spec- 
trum of fluctuations. The best studied example, both 
in theory [2, [3, i, i, H, 0] and in experiment [I,!, is 
the case of two-level fluctuators in the environment of a 
Josephson qubit If the number of fluctuators cou- 
pled to the qubit is small, the decoherence of the qubit 
shows evidence of nonGaussian environmental fluctua- 
tions. Another important class of such environments is 
provided by mcsoscopic transport, in which a quantum 
system is influenced by the nonequillibrium transport of 
electrons through some device. Examples of such en- 
vironments include the partition noise from a quantum 
point contact [l^, j and transport through a single- 
electron transistor (SET) [li]. 

In this article, we describe a general theory of the in- 
fluence of nonequillibrium, nonGaussian environments on 
quantum dynamics. We work within a generalised master 
equation (GME) framework [l^, [l^ , and assume that the 
environmental degrees of freedom to which our quantum 
system directly couples may be described by a Marko- 
vian GME of the Lindblad form. Under this assumption 
we derive an effective Liouvillian describing the reduced 
dynamics of the system alone that can be expressed in 
terms of environmental correlation functions. We thus 
obtain an explicit account of the effects on the system of 
environmental fluctuations of all orders. 

Whilst this theory is presented in general terms, we 
specifically have in mind applications in mcsoscopic 
transport, where we seek to describe the dephasing and 
relaxation of a quantum system due to the charge fluctu- 
ations of a nearby mcsoscopic device. In Ref. [1J| the be- 
haviour of a charge qubit was related to the charge-noise 
spectrum of SET environment. It is one of the aims of 
this work to place the results of Ref. in a broader con- 
text and to generalise not only to arbitrary mcsoscopic 
devices in the Coulomb blockade regime, but also to in- 
corporate the effects of chargc-fiuctuations of orders be- 



yond Gaussian. For a qubit coupled to the environment 
via a pure-dephasing coupling, we describe how the long- 
time behaviour of the qubit is related to the cumulant 
generating function of the operator through which the 
system couples to the environment, and show how to cal- 
culate this quantity for arbitrary environments. 



As illustration of this theory we consider a charge qubit 
couplied to two mcsoscopic environments: i) the SET en- 
vironment of Ref. which is equivalent to a source 
of classical telegraph noise [T7|. and, in certain limits, a 
model of a single background charge fluctuatorQ, 01, and 
ii) a double quantum dot (DQD) environment. Whilst 
transport through a DQD has been extensively studied 
[l6l . [isl . [19I . HO, El, to our knowledge, its role as a deco- 
herence source remains unexplored. Moreover, the DQD 
environment is an important example because, whereas 
the SET model can be described in purely classical terms, 
the inter-dot coherence of the DQD means that trans- 
port through it, and hence the fluctuations to which the 
charge qubit couples, are quantum mechanical in nature. 
Both these examples exhibit interesting nonMarkovian 
qubit dynamics, including dramatic visibility oscillations 
in the strong coupling limit . The DQD model also 
exhibits a pronounced quantum Zeno effect |22| in this 
same limit. Comparison of these models highlights the 
distinctions between quantum and classical fluctuations 
in determining the dephasing and relaxation of a system 
coupled to them. 



This paper is organised as follows. We flrst describe 
the general model considered here and its description in 
terms of coupled GMEs. We then show how the envi- 
ronment may be traced out to a yield an effective Li- 
ouvillian for the system. This Liouvillian is related to 
environmental correlation functions, used to derive de- 
phasing and relaxation rates for the system in the weak 
coupling limit. Two special cases are then discussed in 
which the results are particularly simple: pure dephasing 
and classical environments with relaxation. We conclude 
with a study of our two examples and discussions. 
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FIG. 1: Quantum system S is coupled to an environment 
divided into two parts: E, which couples directly to S, and 
E', which does not. E' is large, in equilibrium, and weakly 
coupled to E such that E is maintained in nonequilibrium. 



I. SYSTEM-ENVIRONMENT MODEL 

Figure [T] depicts the general situation under discus- 
sion here. The environmental degrees of freedom are 
divided into two sets, labelled E and E', according to 
whether they couple to quantum system S directly or 
not. The Hamiltonian of the system-environment com- 
plex isH = Hs+'HE + He' + gVsE + Vee' , with 7Ys,e,e' 
the isolated Hamiltonians of our decomposition, Vse and 
Vee' interaction Hamiltonians between system and envi- 
ronment E, and between environmental components, and 
g a dimensionsless coupling constant. We assume that 
the EE' coupling is weak, that reservoir E' is in equilib- 
rium, and that the Born-Markov approximation is valid 
for the EE' coupling. Following a standard master equa- 
tion derivation we trace out environment E' and obtain 
a GME for the SE density matrix: 



dtp' 



SE 



S^S 



Q + gM 



SE\ ^SE 



(1) 



with system Liouvillian CqP' 



i [Hs,P^] , SE couphng 
LiouviUian A4p^^ = — i [Vse,P^^], and LiouviUian Cq 
given by a Lindblad form obtained by tracing out E'. 

In the following we will employ a notation for GMEs 
in which the elements of the density matrix are arranged 
into a vector \p^^)) with populations first, followed by co- 
herences [2^. In this notation, superoperators are writ- 
ten as matrices, and the GME for the SE density matrix, 
now a vector, has the form 



dt\p' 



SE. ^ /:SE|^SE^^^(^S + ^E+^_^)|^SE^^ 

= {Cl^+9M)\pn). (2) 



The GME for the environment reads dt\p^)) = £^\p^)). 
Let us denote the eigenvalues of Cq as X^, assumed dis- 
tinct, and its right and left eigenvectors as and 
respectively. These vectors form a biorthogonal set. 



{{4>m\4'n)) = Snm, but are not adjoint, since £q is non- 
Hermitian. The stationary state of the environment is 
given by |p^tat)) = l'?^o))' ^he zero-eigenvalue eigenvector 
of C^, i.e. C^\(l)^)) = Xf\(t>f)} = 0. The corresponding 
left-eigenvector ((0o L tias elements 1 at locations corre- 
sponding to populations and is zero otherwise. Similar 
definitions hold for the free system Liouvillian and its 
eigendccomposition; but note that, since Cq contains no 
damping; its nuUspace will be of dimension greater than 
one. 

Finally we assume that the SE interaction has the bi- 
linear form, VsE = ^cre, where cr is a dimensionless sys- 
tem operator and e an environment operator with di- 
mensions of energy. The corresponding Liouvillian is 



SE 



p^^aej 



which we 



obtained from A^p^^ = —i-^ [<^£p 
write in vector notation as 

Mlp'"")) = -^l {OtOt 0-0-) IpSE)). (3) 

The superoperators Of, here represented as matrices, 
can be obtained by considering matrix elements. 

Although we will derive results for arbitrary systems, 
it is often useful to discuss the case when the system 
is a qubit. In its diagonal basis, the qubit Hamiltonian 
is Rs = 5^0-2 and, in the basis (pn, P22, P12, P21), the 
corresponding free Liouvillian is 
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(4) 



The system part of the SE coupling operator g is then a 
traceless Hermitian matrix with elements n • cr where n 
is a unit vector and <t the vector of Pauli matrices. The 
relevant operators in Liouville space are 
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(5) 



with 7i± = Tlx ± 71-y. We will discuss two particular ex- 
amples in the following: 

Pure dephasing: With cr = CTz, the coupling is diagonal 
in the same basis as the free evolution of the system. In 
this case, we have a 'pure dephasing' model with both 
of diagonal: 



O^ 



/ 1 

0-100 
10 

Vo 0-1 



o: 



/ 1 00 

0-100 
00-10 
Vo 1 



.(6) 



Under this coupling, only the off-diagonal elements of 
the qubit density matrix evolve in time, and we write 
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Poi{t) = D{t)poi{0) with the 'degree of coherence' D{t) 
to be determined. 

Orthogonal coupling: The other couphng that we will ex- 
plicitly consider is 'orthogonal coupling', with a = a^. 
In this case we have the off-diagonal super-operator ma- 
trices 
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(7) 



and this coupling will induce relaxation in the system. 

Finally, we will also refer to the situation when the en- 
vironment is classical. In this case, only populations are 
required to describe the state of environment E, and the 
GME determining its behaviour is actually a rate equa- 
tion. It then follows that operators acting on E are diago- 
nal and commute at different times. This in turn implies 
that the ± superoperators are identical: 0~ = 0+. As 
will be made clear below, this situation represents an en- 
vironment which experiences no back-action due to its 
interaction with the system. 



II. EFFECTIVE SYSTEM LIOUVILLIAN 

Having set up our model of coupled master equa- 
tions, we now proceed to derive a description of the sys- 
tem's behaviour in terms of environmental quantities. 
To this end, we derive an effective Liouvillian for the 
system dynamics. Laplace transform of Eq. ([2]) gives 
|pSE(^))) ^ 17SE(^)|^SE(^ ^ 0))) with SE propagator 

il^^(z) — [z ~ . This we expand in orders of g 



nl^iz) + gnl^iz)Mnl^iz) 
+g'nl^{z)Mnl^{z)Mnl''{z) + ...,{&) 

with rio^(^) = [z - C^^] ~\ the free SE propagator. As- 
suming that the environment starts in its steady-state 
pfjjjt, the reduced system propagator D,^{z) is given by 
il^{z) = |r2^^(z)|0Q )), corresponding to a trace over 
the remaining environmental degrees of freedom. With 
the expansion of Eq. ^ , we have 

n^{z) ^ {{^fin'^izMf}) 

- nt{z) + nUz)U^\M\ci)^))nliz) 
+n^oiM^o\Mnl^iz)M\4>^))nUz) + . . .(9) 

where we have identified the free system propagator 

nliz) = {{c^f\n^^(z)\^^)) = ^—. (10) 
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We can consider the system propagator as arising from an 
effective system Liouvillian, C^g{z), which will be non- 
Markovian. We write this as a series in g: 



-oft 



E 

n=0 



(11) 



such that the full system propagator can also be written 
as 

c 1 

n^{z) = 



- - ^fft 

= VlI{z) + n%{z)L\VLl{z) + VlI{z)11%VlI{z) 
-rVLl{z)Z\VLl{z)L\VLl{z)^.... (12) 

An order-by-order comparison of Eq. ^ and Eq. (flT 
gives us 

Formal resummation yields 



' M\ 



n>\. (13) 



1 



A^I</'o))Xl4) 



where 1^ — I0o ))((<^o L ll^*^ projector out of the en- 

vironment steady state. Using the eigendecomposition of 
the free SE Liouvillian, we can write the free SE propa- 
gator as 



^l^{z) 



1 



E 



1 



E ; - AS - AE 



(15) 



with Zu — z — The effective Liouvillian may then be 
written 



■^off 



(^) = 
+5((0o 



-^0 



(16) 



The effective Liouvillian of Eq. IH]) or Eq. HI]) is the 
main formal results of this work; it describes the system 
dynamics in a compact, self-contained form and includes 
environmental fluctuations of all orders. In this form it 
is not particularly instructive, however, since both S and 
E quantities appear in an intertwined way. In order to 
see the significance of these results then, we will consider 
first a weak coupling expansion, and then some special 
cases where S and E dependencies can be separated. 

Before doing so, let us note that an expression similar 
to Eq. ^ can be written down for the reduced envi- 
ronmental propagator, f2^(z). The initial state of our 
system, which we take to be a qubit here, is arbitrary 
and can be written in the form 



\pI)) 




(17) 
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which defines the vectors |0f )), and coefficients — 1 < a < 
1 and b < 1 — a^. The only conjugate state we shall need 
here is ((^qI (1,1,0,0). With the system starting in 
this state, the environmental propagator can be written 
as 



= n^{z) + gn^{z)U^,\Mnl^{z)\pl)) 



(18) 



with Se(t) = e(t) — (e). Here the time-dependence of the 
operators is given by the evolution of the full environmen- 
tal Hamiltonian He + We' + Vee' • Using the QRT to ex- 
press this correlation function in terms of quantities act- 
ing on E alone, we obtain S^^) (z^) = {(0+n^iz^)Q^O+)). 
Similarly, by recalling that super-operator 0~ is equiva- 
lent to operator e acting from the right, we obtain 



{{otnf{z,)Q^ot)) = s^'Hz^) 
{{otn^iz,)Q^o-)) = (5(2)(z:) 



(21) 



With this expression we can calculate the effects of back- 
action of the system on the environment. We will not 
follow this calculation further, except to note what hap- 
pens for classical environments. In this case, we have 
O- = 0+ and thus M = iO+ (0+ - O"). The ef- 
fective environment propagator thus contains terms like 
iO+{{(j)^\ (0+ - 0-)n^'^{z)\p^)) and, as is easy to ver- 
ify, ((0§|(O+-O^) = 0. AU terms beyond the first 
in Eq. p8)) start in just this fashion, and therefore, 
for classical environments we have 51^ (z) = fig (z), and 
there is no back-action. It is possible to construct 
models with a classical environment that do experience 
back-action with, for example, the rates of the free en- 
vironmental Liouvillian depending on the state of the 
qubit [111, [13, t2^] Such models, however, are outside 
the class discussed here in which all back-action effects 
arise from the quantum-mechanical nature of the system- 
environment coupling. 



III. WEAK COUPLING: DEPHASING AND 
RELAXATION RATES 

We now consider the situation where the SE coupling 
is small g ^ I and describe the weak coupling expan- 
sion of Eq. (jl6p . In this case, we can consider the par- 
tial LiouviUians of Eq. (|13|) as successive approxima- 
tions to the full Liouvillian. At first order, we have 
Cf = i(e) (0+ - 0-), with (6) = ((0+)) = ((Or)) the 
steady-state expectation value of operator e. 

At second order we have 



'-2 



U^\MQ''nl''{z)Q^M\<^^)) 
E«'^ol>'l€))((^^|Q''^^o(^.)Q'=-^l0o))-(19) 

v=0 



From the form oi M, we see that this expression depends 
on environmental quantities like ^O'^ Q^^l^{z^)Q^Of)) , 
which can be evaluated straightforwardly for any partic- 
ular model. Moreover, they can be related to correla- 
tion functions of operator e via the quantum regression 
theorem (QRT) [25l. [26j. Let us define the second-order 
correlation function 



dre 



'{Se{r)6ei0)) 



(20) 



Putting these results together, we obtain our final form 
for the second-order effective Liouvillian 



{~^g/2rJ2iOt-0-)\c|^l))m 

1^=0 

x{0+5(^)(z.)-0,"(5(^)(z:))*}. (22) 



This Liouvillian determines the system behaviour for all 
times in the weakly coupled limit for arbitrary system 
and environment. Its form is simply that of a matrix in 
system-space, the elements of which contain environmen- 
tal correlation functions evaluated at various frequencies. 

The long-time behaviour of a qubit can be described 
by a pair of rates, Td and describing dephasing and 
relaxation respectively. These rates are determined from 
C^ff{z), or as is the case here, its second-order approxi- 
mation. We first diagonalise C^ff{z) — V{z)A{z)V~^{z), 
such that the effective system propagator may be writ- 
ten = V{z)j^j^V-^{z). We then find the poles of 

r, with A;;fc(z) the kth eigenvalue of C^ff{z). 



each — r — 7 — r . 

z-Akk (z) • 

In the long time limit, only the pole lying rightmost in 
the complex plane contributes. For a qubit then, ^ 
assumes the asymptotic form 
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Vo 
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C2 







J 



(23) 



with Ci constants, ly some frequency of coherent oscilla- 
tion and Td and F^, the aforementioned rates. 

We now discuss Eq. and the corresponding rates 
for two illustrative couplings. 

Pure dephasing: Let us consider a qubit and set (e) = 0, 
since we can always incorporate it into a redefinition of 
A. The effective Liouvillian of Eq. (p2)) is diagonal, has 
zeroes at first and second diagonal elements, and has the 
third element 



= l{z) = ^A - {S{z - ^A) + S*i-z + lA)} 



and (£off)44 = (£off)33- The off-diagonal elements 
of the qubit density matrix there evolve as poi{z) = 



5 



D(z)poi(i = 0) with D{z) = {z - l{z)} \ The 
rightmost-lymg pole of D{z) is 

zo = tA-^g^SiO) + S*{0)} , (25) 

correct to second order in g. The dephasing rate is there- 
fore 

Td^lg'si'Ho), (26) 
with the fuU, symmetrised correlation operator 

= 2 J ^dre^^^{{Se{r),Sem), (27) 

with {...,...} denoting the anticommutator. 
Orthogonal coupling: We now consider the orthogonal 
coupling and set (e) ~ for simplicity's sake. The non- 
zero poles of the second-order effective Liouvillian with 
system operators of Eq. ^ are 

zi = -(g/2)2{5(-^A) + 5(^A) + 5*(-^A) + 5*(^A)} 

Z2 = iA-{g/2f{S{iA) + S*i-iA)} 

Z3 = -lA- (5/2)2 {5(-^A) + 5*(^A)}. (28) 

The corresponding relaxation and dephasing rates are 
r.r = zi = ^g^S^iA); = Re (23) = ^T,, (29) 

correct to second order in g. The stationary state of the 
qubit is found to be 

/ SiiA)+S'iiA) „ \ 

/ - - - (30) 

Pstut — I S{-iA)+S'{-iA) ■ y'^^) 

V " sf'(A) J 

We define 

P = 2Tr{y92}-l (31) 



as a measure of the purity of a qubit density matrix. 
For a completely mixed state P = and pure state 
P = 1. For a pure dephasing model, the final purity 
depends on the initial state — the generic final state is 
p = Diag (c, 1 — c), where c is given by the initial condi- 
tions. The purity of this state is P = (2c — 1) , which is 
unity if we start in a pure localised state, c = 0, 1, and 
zero if we start in the superposition c = 1/2. For the sta- 
tionary state of the orthogonal coupling model, Eq. (|30p . 
the purity is 



^ = ( TPTTTT r dre'-^iMr), <5e(0)]) ) (32) 

\Se (A) J-oa j 

with [...,...], the commutator. In words: the pu- 
rity is determined by the ratio of the Fourier trans- 
forms of the commutator and the anticommutator of 
fiuctuation-operator (5e at different times. For a classi- 
cal environment, the commutator is zero, and the pu- 
rity is zero. Nonzero values of the purity are an indica- 
tor of the existence of back-action of the system on the 
environment [31 , and in the current class of models, the 
back-action is always quantum. 

These results illustrate the generality of the connex- 
ion between these rates and the second-order correlation 
functions given in Ref. [l3 |. In principle, we can extend 
the above analysis to arbitrary order in making the 
connexion to environmental correlations functions with 
the QRT. At third order, for example, we have the effec- 
tive Liouvillian 

L% = Ul\MQ''^l^[z)Q^MQ''^l^[z)Q''M\<^l)). 

(33) 

With third-order correlation functions defined as in ap- 
pendix]^ this Liouvillian can be written as 



-o~m)ui\ (ot (s(^'\z:„zi)y -o- (5f")(z:,z:o)*)}- (34) 



IV. PURE DEPHASING 

The Liouvillian of Eq. contains a mixture of sys- 
tem and environment operators in the inverse and in gen- 
eral this means that this inverse can not be carried out 
explicitly. Useful results can be obtained by expansion, 
as above, but the results can becomes unwieldy. In cer- 



tain cases, however, we can effect a separation of S and 
E components, expressing the matrix elements of C^q in 
terms of environmental quantities. In this section and 
the next, we consider two cases: pure dephasing and or- 
thogonal coupling to a classical environment. 

In the pure dephasing model, the effective Liouvillian 
is diagonal, and only has non-zero elements at positions 3 
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and 4. The degree of coherence D{z) is given by the third 
diagonal element of fl^ (z) . In an interaction picture for 
the system (obtained by shifting z — » z — £q), we obtain 



D{z) 



liz)Y 



(35) 



with 



l{z) = -ig{{<p, 



igO,n''{z)Q] 'O,|0^)), (36) 



where we have used the shorthand Oe ~ (O^ + 0~ ) . 
This quantity contains information on environmental 
fluctuations of all orders, as may be seen by expansion. 
However, expansion is not necessary because, for a finite- 
dimensional environment E at least, it may be evaluated 
directly in closed form through matrix inversion. Eval- 
uation for an infinite-dimensional environment may still 
be possible using phase-space methods p6j . 

Equation Eq. ([55]) also suggests the point-of-vicw of the 
purc-dephasing qubit acting as a detector of the environ- 
mental fluctuations. As we now show, D{t), the inverse 
Laplace transform of D{z), is related to the generating 
function for the zero-frequency cumulants of operator e. 
The inverse Laplace transform of D{z) is obtained from 
the roots of z — ^z) = [l^]. In the long time limit, 
we only need the pole lying rightmost in the complex 
plane, which we denote z, such that we can approximate 
D{t) = ce^* with c some constant. Considering g as in- 
finitcsimally small, we can write z as an expansion about 
z = 0: z = X]n=i(~*5)"'^n ^^'^ solve z — l{z) = order- 
by-order. For example, the first three terms are 

Zi = ((O,)), f2 = ((af^^(o)QOe)), 

Z3 = ((O,f^^(0)Q(O, -((a)»^^^(0)QOe)), (37) 

with the shorthand ((...)) = {{<Po\ ■ • ■ |0o ))• Application 
of the QRT, shows that these quantities are equal to the 
zero-frequency limits of the Keldysh-ordered correlation 
functions: zi = (e), Z2 is equal to the zero- frequency limit 

of Sli^^uj) from Eq. ([TT)) . and Z3 is the zero- frequency 
limit of S'P' (w, w') defined in appendix [XI and so on. 
Therefore, within the approximations made here, z may 
be identified as the cumulant generating function for 
zero-frequency, Keldysh-ordered correlation functions of 
the operator, which we denote Se, Se"\{0}). The cou- 
pling parameter g is identified with the 'counting field' 
such that 



Qr, 



d{~igy 



5(")({0}). 



(38) 



9=0 



This result mirrors that of full counting statistics 
[Hi [13 in which the dephasing of a probe qubit is re- 
lated to the cumulants of the number of electrons passed 
through the device, viz. the zero-frequency correlation 
functions of current fiuctuations 51. In this context, the 
relationship between Keldysh-ordered CGF and dephas- 
ing was to be expected. However, by proceeding in the 
above manner, we have obtained explicit expression for 



the full CGF, as the inverse Laplace transform of D{z) 
with l{z) as in Eq. (|36p . as well as for the lowest cumu- 
lants from Eq. ([57]) . Obtaining these expressions from 
standard Keldysh approach would be non-trivial. The 
simplicity arises here from the introduction of the two su- 
peroperators O^, which correspond to the two branches 
of the Keldysh contour. 

We have talked considerably about the long-time limit 
of the behaviour of the qubit. However, the short-time 
behaviour is also of interest, particularly with an eye to 
applications in quantum information processing where 
we definitely want to avoid the heavy dephasing of the 
long-time limit. In such circumstances, only the weak- 
coupling limit is of interest and thus we restrict our short- 
time discussion to the second-order approximation. We 
have D{z) as before, but from truncating Eq. p6p . we 
approximate 



l{z)^^g^{^^\0,n^{z)Q^OM^)), 



(39) 



with first-order contribution removed by considering the 
appropriate rotating frame. As shown in appendix [Bl we 
can explicitly perform the Laplace transform of the D[z) 
arising from Eq. (|39p and, correct to order g"^ and leading 
term in i, the coherence decays as 



(40) 



The leading term is proportional to , which is to be 
compared with the linear dependence that would arise 
from a simple Markovian decay D{t) ~ 1 — F^f. Note 

also that it is the quantity \g'^{{(j)^\0^ Q^Od^Q )), and 
not Fd, that determines the time-scale of this initial be- 
haviour. 



V. ORTHOGONAL COUPLING AND 
CLASSICAL ENVIRONMENT 

We can also effect the separation between system and 
environment quantities if environment is classical. For 
simplicity, we consider just the orthogonal coupling of 
Eq. ([7]). To proceed, we define the system operator 



F 




(41) 



and use it to transform the effective system Liouvillian 
of Eq. mi). With C'^g = F£fffF-\ we find 



f - 
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A2 















Vo 




•^44 / 



(42) 
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with elements 

^24 = ^9{{<P^\AOtm 
= ig{{4>^\B0t\^^)) 



-42 

^4 = -5'((</'^|i3O+17oSO+|0^)), (43) 



operators 



and 



17j 



1 



{l7[;^(z + iA)±17^(z-iA)} 



- l^S'l-^)- (45) 

The stationary state in the original basis is Pstat = 
i (|0}(0| + |1)(1|), the fully incoherent superposition of 
the two qubit states. This holds true for any non-pure- 
dephasing coupling to any classical environment without 
back-action. 

VI. SINGLE ELECTRON TRANSISTOR 

As a first example, we consider a charge qubit coupled 
to a SET environment as discussed in Ref. The en- 
vironment consists of a single level coupled to electron 
reservoirs that can either be occupied or empty. The 
infinite bias limit is assumed such that transport is uni- 
directional, with electrons entering the level with rate 
and leaving with rate F/j. The only pertinent density 
matrix elements are the two populations of an empty or 
occupied level. In this empty/full basis, the Liouvillian 
of single level reads 



- 



(46) 



The stationary state of the environment is |0o)) ~ 
{VR,TLf , with r = Tl + Vr, and £f has a sin- 
gle non-zero eigenvalue: Ai = —T. The charge qubit is 
taken to couple to the charge in the SET, and we take 
the coupling operator to be e = FrtgET with tt-set, the 
number operator for the single-level, and F setting the 
energy scale. The two coupling superoperators are 



ot = o: 



F 




1 



(47) 



which are observed to be equal since we have a classical 
environment. 

With the qubit coupled to this environment in a pure 
dephasing configuration, the effective Liouvillian term 
l{z) of Eq. pS]) evaluates 



l{z) 



-«5 



Fl (F + z) 
{T + z) + igTi 



(48) 



The simplicity of this model means that the inverse 
Laplace transform of D{z) = [z — l[zy)~^ can be per- 
formed analytically. With F^ — Tr for simplicity, we 
have 



{cosh(CFfli) + C"' sinh(CFflO} 



(49) 



with C, = \Jl — . For weak coupling, g <C 1, this ex- 
pression correct to second order, becomes 



(44) D{t) w e-'f rH*e-^* ( 1 + ^ ( 

2Fd V 



-2(rH-7)n 



)|,(50) 



with weak-coupling dephasing rate 



In the strong coupling limit, g 3> 1, we obtain 

D{t) « e-'f^«*e-^«* {co&{vt) + g-^ &in{vt)] 
with oscillations at a frequency 
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(51) 



(52) 



(53) 



and a decaying envelope of rate F^ = F^. These two 
contrasting behaviours are illustrated in Fig. [2^ in which 
we plot the visibihtyfH] v{t) = \D{t)\. Our approach 
here therefore reproduces the salient features of the pure- 
dephasing results for a classical telegraph noise model [3, 

Taking the long-time limit of Eq. , we obtain the 
CGF for the SET charge operator 



rj? (-i.g-l + VT" 



(54) 



The first few correlation functions are determined from 
Eq. ^ as Tr (1, 0, -3, 0, 45, ... - 1575, 0, . . .) such that 
the odd-functions above the first are zero, and are seen 
to grow with increasing order fc as ~ . 

With the SET coupled to the qubit with a cr^ coupling, 
S and E contributions can be separated along the lines 
of section |Vl With F^ = Tr, the effective Liouvillian in 
the basis of Eq. (j42|) is 



^cff 



/O 

k{z + 2Tr) igVR{\ + K) 
-iA 

A^ + (z+2rH) 



with 



\0 igTR{l + n) iA (j+ar^ 
A2 + (.gFfl)2 + {z + 2Ffl)2 ■ 



(55) 



(56) 



With off-diagonal elements in the effective Liouvillian, 
the behaviour of the qubit is more complex than in the 
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FIG. 2: Left panel: The qubit visibility v{t) as a function of 
time for a SET environment with a pure-dephasing coupling. 
As g increases through g = 1, the behaviour shifts from a 
simple decay to a highly nonMarkovian oscillating one. Right 
panel: Magnitude of the coherence |poi(i)| of the qubit cou- 
plied to the SET environment with the orthogonal coupling, 
starting from an initial state 2~^^'^(\0) + |1)). Here, the weak- 
coupling case g <^ 1 produces coherence oscillations, whereas 
strong coupling induces a simple decay. In all cases, we set 
Tl = Fr, and A = O.STr. 



pure-dephasing case. At first order in g, the environ- 
ment induces a tr^ term in the system LiouviUian. In 
the weakly coupled case, 5 <C 1 , this leads to coher- 
ent oscillations between the qubit states, as Figure [2b 
shows. These oscillations are damped with rates = 
25^ry (A^ + 4F|j) and F^ = Fr/2. As g increases, the 
damping becomes stronger until the oscillations are sup- 
pressed for (7^1, A/Tji. This behaviour with increasing 
g is exactly the opposite as found for the pure-dephasing 
model. 



VII. DOUBLE QUANTUM DOT 
ENVIRONMENT 



We now consider a charge qubit coupled to the elec- 
tronic transport through a DQD [l^, [l^, [l^, [2^ , as de- 
picted in Fig. [3l The crucial difference between this and 
the previous SET environment is that here the position of 
the electron in the DQD is a quantum-mechanical vari- 
able [21I, [s^l, and thus requires a full GME treatment, 
and not just a rate equation. 

We assume the strong Coulomb blockade limit such 
that the only relevant states are the empty DQD, and 
a single electron in either left or right dot. In pseu- 
dospin language, the Hamiltonian of an electron in the 
DQD reads He = ^scrf + Tca^, with detuning e and 
tunnelling element T^. The two levels are coupled to 
their respective leads with rates F^ and F/j, . The Liou- 
viUian for such a DQD[23| may be written in the basis 



FIG. 3: Our second example consists of a charge qubit system 
(labelled and 1), coupled to environment consisting of a 
double quantum dot E (labelled L and R), and environment 
E' consists of two electron reservoirs (labelled and fin). 
The DQD is biased such that it is in the nonlinear transport 
regime. 



{pm, Pll, Prr, Plr, Prl) as 



^0 



-Fl 
Fl 









Ffl 
iTc 

-Tr -iT, 

iTc -iTc -ie - ^Tr 

-iTc iT, 





-iTc 
iTc 


-h^RJ 



(57) 



The coherent level splitting is Adqd = \/e^ + 4T^ and 
in the following we will assume the rates to be equal F^ = 
Tr. We will also assume that the qubit can be coupled 
to the DQD environment through any of the three Pauli 
matrices, as well as through N, the total occupation of 
the DQD. A related model was discussed in Ref. [3l| but 
with transport through both sets of quantum dots, i.e. 
both through the DQD of the environment and the DQD 
of our charge qubit. 

Assuming that we can couple to any of these opera- 
tors, the qubit in pure-dephasing configuration allows us 
to perform a complete tomography of the DQD fluctu- 
ations through an application of the results of section 
IIVI By this it is meant that, since the density matrix 
is completely determined by measurement of the ex- 
pectation values of these four operators, the long-time 
fluctuation spectrum is also completely mapped by the 
correlation functions of the same four operators. In Fig- 
urc[31 we plot this set of correlation functions up to fourth 
order as determined from the generating function z. All 
correlation functions are even functions of e with the ex- 
ception of the odd-numbered functions for cr^. . Operator 
ax also stands out as having fluctuations which grow ex- 
tremely rapidly with increasing order. In the long time 
limit, the expectation value of ay gives the stationary 
current through the DQD, as expected [ll]. Surprisingly, 
although ay can not in general be identified as the cur- 
rent operator through the DQD, the second correlation 
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FIG. 4: The first four zero-frequency correlation functions 
si"'' (0), n — 1, ... ,4, of operators e = Ox, cry, a^, and A'' for 
the DQD environment as a function of DQD detuning e/Tr. 
The tunnel coupling was Tc — Fr, with Fl = F_r. 



function of g,, here does correctly reproduce the DQD 
shotnoise [li,l2l|. 

We now turn our attention to the dynamics of the 
qubit itself, and focus first on the pure dephasing case. 
The first coupling operator that we consider is Vse = 
iPijcrfcr^. In the weak coupling limit, g 1, the de- 
phasing is well described by the expansion of Eq. ([5^ 
up to second order. Figure [5] shows the behaviour of 
the qubit visibility over both long and short time scales. 
The dephasing rate F^ ~ Ixi^) is easily calculated from 
Eq. (|26p . but its form is too unwieldy to give here. In the 
special case of e = 0, however, it reads — ^ 0) ^ 
2g^rR{Tl + 8r2)/(F^ + UT^/Tr). The Corresponding 



visibility as a function of time is 



27. 



-(rH/2-7x)t 



(58) 



which, at short times, reads v{t) = 1 — jx^Rt'^ /2, show- 
ing the generic deviation from Markovian behaviour. For 
finite e, the expression for the visibility is more compli- 
cated and, as Fig.[5]shows, describes oscillations at short 
times. This nonMarkovian effect is due to the coherent 
coupling of the system to the environment, with the oscil- 
lations proceeding with a frequency ^ Adqd- It should 
be noted that these oscillations will only be visible if the 
DQD dynamics are not too heavily damped pol |. 

The e = case is a particularly special limit of this 
model in which the environment is, at second-order in 
g, identical to the SET, as a comparison of Eq. ([58]) 
and Eq. (|50p shows. The reason for this equivalence 
is most easily seen by rewriting Eq. ([FT)) in the basis 
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FIG. 5: The visibility v{t) of the charge qubit coupled to 
the DQD environment with a^a^ coupling. Lefthand panels: 
Weak coupling behaviour with g — 0.2 for different values 
of e/Tr. The time axes are scaled with weak coupling rate 
^x{e). On longer timescales (a), the behaviour is approxi- 
mately Markovian but at short times, deviations are apparent 
(b). Righthand panels: The behaviour for strong coupling is 
highly nonMarkovian, with pronounced visibility oscillations. 
SE coupling strength and DQD detuning were here g = 8, 
£ = (c) and g = 4, e = 2 (d). The time axes are scaled with 
Tr. For all plots, we set Tc — Tr. 



(poo, Pll,Pbm, Re(pLfl), Im(pLfl)): 
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It is then clear that for e = 0, the real part of the inter- 
dot coherence, Re(pLij), decouples from the rest of the 
(free) environment dynamics. At second order, the cou- 
pling operator a^'^^ couples only to this real part of the 
coherence and the DQD environment therefore appears 
as a classical environment with a single excited level of 
width Tr/2. At orders higher than second, this equiva- 
lence breaks down as the operator couples to further 
parts of the environment density matrix. 

For stronger coupling g ^ 1, nonMarkovian effects 
dominate. We consider first the case with e = 0. For 
g ^ 1 the visibility behaves like that of the strongly- 
coupled SET with oscillations at frequency ^ g^R and 
an envelope that decays as Tr/2. Increasing g further 
leads to a doubling in the period of the oscillations not 
seen in the SET model. This behaviour can be traced to 
the increased importance of a further pair of poles of D{z) 
oscillating at a frequency ^ 2Tc for g ^ 1 . As illustrated 
by Fig. [SJi, the situation for e 7^ is more complicated, 
with frequencies associated with all coherent couplings of 
the SE complex playing a role. 

Now consider a second SE coupling operator: Vse = 
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FIG. 6: The visibility v{t) as a function of time for the cou- 
pling VsE = ^r_Rcrf (T^. The main panel shows that the visi- 
bility decays with an initial rate rii/2 until it reaches a value 
which appears stationary on the scale of this plot. Also shown 
is the population of the left dot of the DQD, pll, which shows 
localisation due to the Zeno effect induced by the coupling to 
the qubit. The inset shows the visibility on a much longer 
time-scale, which reveals a slow decay with rate 7zcno. The 
parameters were g — 12, Tc — Tr and e = 0. 



^Tjiafaf. For weak coupling, the behaviour is qualita- 
tively similar to that for the previous CT;^-coupling, with 
the £ = dephasing rate given by 72 = 4g2r|.T^^(13r|.+ 
36T2)/(r| + 12T^f. The strong coupling behaviour 
is, however, significantly different, as can be seen from 
Fig. [6l Initially, the visibility shows oscillations and de- 
cays with rate Tji/2. But instead of the visibility de- 
caying to zero after a few times F]^^, it reaches a value 
which, on the time scale of the initial decay, appears con- 
stant. This is a manifestation of quantum Zeno effect, 
with the coupling to effectively measuring and local- 
ising an environmental electron in the Icfthand dot of the 
DQD. Transport through the DQD, and along with it the 
dephasing of the qubit, is then almost completely inhib- 
ited. For large but finite coupling, this localisation is 
not perfect, but decays with the rate 7zono = T^/ig^^n) 
for e = 0. Needless to say, this effect relies on having a 
quantum environment. 

Finally, we consider couplings to the qubit such that it 
undergoes not only dephasing but also relaxation. The 
behaviour is correspondingly more complex and here we 
focus on just one aspect: the effects of quantum back- 
action. In section |V] we saw that for a qubit cou- 
pled with any non-pure-dephasing coupling to a clas- 
sical environment, the stationary qubit state is always 
Pstat ^ ? + Let us therefore consider a 

SE coupling VsE = ^^R'^x'^x ^"^^^ focus on the weak- 
coupling limit. Fig. [7] shows the behaviour of the popu- 
lations of the two qubit levels as a function of time for 
e = Ffl/2. As is clear, a steady state is reached in which 
Poo 7^ Pill or in other words for which (tr^) is in the 
steady state is non-zero. Due to the dephasing of off- 
diagonal elements, the purity of the final state is sinply 




FIG. 7: DQD model with crfcr^ coupling. Left; Popula- 
tion of the two qubit states as a function of time starting 
from an equal superposition with e = Tr/2, g = 0.4 and 
A = Tc = r_R. A population difference develops due to the 
quantum nature of the environment. Right: The stationary 
qubit population difference (irf ) as a function of e for various 
values of Tc/Tr. For e = 0, we have (erf) = and the re- 
laxation of the qubit is effectively classical. Away from this 
point, quantum back-action leads to a finite value of (cr?). 
Parameters were g — 0.2, and A = Fij. 



given by P = {az)'^. In Fig. [TJd we plot (cr^) as a func- 
tion of detuning, from which it is immediately clear that 
(az) = only if e = and that (ctz) is antisymmetric as 
a functions of e. 

This behaviour is evidence of the back-action of the 
system on the environment. For e = 0, we effectively have 
a classical model in the weak-coupling limit, as discussed 
above, and thus (ct^) = follows accordingly. Even for 
large couplings, this model always relaxes to the com- 
pletely mixed state for e = 0, implying that all coupling 
between system and environment are essentially classi- 
cal and back-action free. For e 0, the coupling to the 
qubit disturbs the motion of the environment which in 
turn feeds back into the behaviour of the qubit and leads 
to a finite value of (ct^) in the stationary limit. 

The relevance of the foregoing results for experimental 
quantum dot systems can be assessed as follows. The 
weak-coupling Markovian dephasing/relaxation rates for 
the models considered here are of the order g^T, with 
r some typical environment rate and where, in princi- 
ple, g can be as small as one wishes. Finite detector 
bandwidth should therefore not be a problem, provided 
sufficiently a weak SE coupling can be realised. Con- 
versely, the nonMarkovian effects occur on energy scales 
of either F or at the Bohr frequencies of the environment 
A^. In early DQD experiments, such as Ref. [s^l, the 
tunnelling rates and coherent coupling strength were of 
the order of ~10MHz. Observing the nonMarkovian ef- 
fects produced by an environment with such time-scales 
would be very challenging. However, more recent exper- 
iments [3^ . [33! operate with coupling parameters in the 
range 1-lOkHz. At such frequencies, the passage of a 
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single electrons through a quantum dot can be detected. 
In a set-up with such parameter, we would expect the 
nonMarkovian effects described here to lie within the ex- 
perimentally observable range. 

VIII. SUMMARY 

In summary, we have presented a formalism for inves- 
tigating the dynamics of quantum systems coupled to 
noncquillibrium environments that exhibit a full range 
of fluctuation statistics. The end product of this analy- 
sis is an effective system Liouvillian which describes the 
systems behaviour without need to explicitly follow the 
environmental dynamics. 

With a weak-coupling expansion, this effective Liou- 
villian can be expressed in terms of environmental corre- 
lation functions, showing how nonGaussian fluctuations 
impact the system behaviour. This represents a system- 
atisation and generalisation of calculations of Ref. [l3| 
to arbitrary systems and environments, as well as to en- 
vironmental fluctuations beyond Gaussian. Away from 
the weak coupling limit, two situations distinguish them- 
selves where the dependence of the systems behaviour on 
environmental correlation functions of all orders can be 
made explicit — pure dephasing coupling and when the 
environment is classical. The pure dephasing case is in- 
teresting as here the qubit can be thought of as a detector 
of the environment fluctuations — a measurement of the 



dephasing of the qubit delivers a generating function for 
the statistics of environmental operator through with the 
system and environment are coupled. 

We have considered two example environments from 
mesoscopic transport — the SET and DQD, acting on 
a charge qubit. These examples serve to illustrate the 
differences between quantum and classical environments, 
the latter class expanded to include situations, such as 
we saw with the zero-detuning DQD with (T^-coupling, 
where a quantum environment acts as if it were classical. 
With a pure dephasing coupling, weak coupling results al- 
ways show deviations from Markovian behaviour at short 
times, but only the quantum environment has the ca- 
pacity to induce oscillations in the visibility. At strong 
couplings, the behaviour of qubit highly non-Mar kovian, 
with the behaviour in a classical environment quite sim- 
ple, but that in a quantum environment far more complex 
with phenomena such as the quantum Zcno effect occur- 
ring. Finally, with a coupling that induces relaxation, 
quantum environments distinguish themselves from clas- 
sical ones, by inducing a steady state of the qubit system 
with a non-zero purity. This is attributed to the quantum 
back-action of the system of the environment. 



APPENDIX A: THIRD-ORDER CORRELATION 
FUNCTIONS 



The third-order contribution to the effective system Liouvillian is expressed in Eq. ((34|) in terms of the following 
two correlation functions 

= ((0+l)^(z)QO+l^^(z')QO+)) 

= / d{ts - h) / d{t2 - ti)e-^(*3-*2)-^'(*2-ti) {{Se{t3)Se{t2)5e{h)) - (e) {6e{ts)5e{h))} 

JQ JQ 

Sf''\z,z') - {{OtVi''{z)QOtVi''{z')QO:)) 

d{h - t2) / d{h - ti)e--(*3-t2)-.'fe-ti) {{5e{t,)5e{H)5€{h)) - (e)(<5e(ti)5e(i3))} 
Jo 

(Al) 

In the long time limit, the generating function z of section HVl is comprised of Kcldysh-ordcrcd correlation functions. 
The third-order correlator [3^ is 

/OC />oo 
d{h~h) d{h~h)e^''^''-'^^+'^'^''-''^s'i\hMM)- (A2) 

with 

S'i\hM,h) - r^(((5e+(t3) + ^e-(t3))('5e+(i2) + <^e-(t2))(<5e+(ii)+<5e_(ti))), (A3) 

where Tx. is the Keldysh-ordering symbol and subscripts ± correspond to forward and backward branches of the 
Keldysh contour. 

APPENDIX B: SHORT TIME, WEAK COUPLING with 

, , . , , , c, = u^\oMl))moM^))- (B2) 

With a pure-dephasmg coupling, the coherence factor 
correct to second order in g is 

1 
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and Oe = ^ (O^ + ). Introducing = g'^Ck/Xf. and correct to second order in g and leading term in t. 
writing the dephasing rate of Eq. (|26p explicitly as = 
-g^Hk^QCk/y^k = -J2k^o'^kXk, D{z) can be written 
in a partial fraction decomposition as 



D{z) 



1 + E 



Oik 



ak 



^ + f-i z - + ctk 

k^O 



(B3) 



It is then trivial to perforin the inverse Laplace transform 
to obtain 



Y^akC 



(Afc-Qfc)t 
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